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Abstract 

Polaward heat transport by the ocean inplles warn water nass fotnatlon, 
l.e. the retention by the tropical and subtropical ocean of sene of its net 
radiant heat gain. Under what conditions net heat retention besostes conparable 
to latent heat transfer to the atnospher^ depends on the relative efficiency 
of transfer processes across the air-sea interface, the top of the atmospheric 
mixed layer, and the floor of the oceanic mixed layer. A thermodynamic model 
of the interacting atmospheric and oceanic mixed layers, with the top of the 
atmospheric layer tak<%n to be at cloud base, shows that net oceanic heat 
retention is significant under the following circumstances. 

-2 

(1) Seasonal heat storage, an^litude of order 100 w m . This is a 
fairly straightforward consequence of the large heat capacity of the oceanic 
mixed layer and leads the seasonal forcing by about a month. 

(2) Massive upwelling (vertical velocity of order 10 ^ m s S , mostly 
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along equatorial cool tongues, with net heat retention of order 100 W m 
The upwelling cooler water is heated and tremsported &vay mainly by the diver- 
gence of surface layer flow (less by the Increasing temperature in the direction 
of the flow) . 

(3) Cold water advection, mostly within the subtropical gyres, net heat 

_2 

retention of order 30 W m . The latitudinal veuriation of radrant heating, 
and generally equatorward surface flow in the northern portions of subtropical 
gyres leads to a moderate rate of warming of the water column as it moves 
along, i.e. to net heat retention of the above order. 

A comparison of model results with observation shows that, over the sub- 
tropical gyres, observed temperature and humidity relationships can only be 


simulated rMliitieally if cold wat«r «dv«ctlon it ttktn into Account. In 
Addition, it it ntcAttAry to tuppott thAt the trAntfor coofficiont At the top 
of thA Atnotphtric mixed leytr (At cloud bASt) it About a:i lergo At At the 
AAA turfACA, whilA thA trtntfAr coAfficitnt At the ocAAnic mixed lAyer floor 
it ntgligiblA, except in regions of mAstive upwelling. The goner a 1 dooiinAnco 
of lAtent heet trenefer erieec from the lerge vaIua of a nondimentionAl lAtent 
heet coefficient (a meteriAl property) end from the repid drop of tAturetion 
cpecific huirddity with height in the Atmosphere. 
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Introdttctlon 

It is now widsly sec^tsa that ths oessn transports hast polswsrd at a 
rata that at cartain latltudaa it of tha ordar of 1 patawatt (10^^ H) , taa 
a. 9 . Hastmrath (1982). Tha plobal proeaaa ia illuatratad acdMOWtically in 
Pig. li in ordar to balanca tha nat haat loaa of tha polar oeaan» watar at 
an avaraga taaiparatura ^ ^ tranafarrad polaward (aithar by "Man" 

flow or by "addiaa") » at tha rata M(kg a*^] . Proai tha availabla taoparatura 
and salinity diffaraneaa and haat and fraahwatar transports it My ba shown 

8 3 ml 

(StORMl and CsAiiady, I960) that N/p ia of ordar 10 a a . h similar 
raault is obtainad by diract oeaanic haat transport and haat budgat calculations 
(Brydan and Hall. 1980} Wyrtki. 1961). 

At high latitudaa. tharafora. %iatar is coolad at tha aggragata rata of 
N kg a , a procass that rasults in "watar Msa formation", i.a. a major changa 
of tha tamparatura and salinity charactaristics of tha watar massas involvad. 

Tha machanism of this procass is fairly wall undarstoodi it involvas vigorous 
panatrativa convaction, vartical mixing, fraquantly. but not always, a daacant 
of tha nawly foxMd mixtura to som subaxirfaca laval and its avantual ascapa 
aquatorward. 

A countar~procasa of warm watar mss foTMtion nacasaarily axists in tha 
tropical ocaan in ordar to raganarata watar massas of highar tamparatura and 
salinity than antar from highar latitudaa. This procass is lass Mil undar- 
stood. Xn contrast to surfaca cooling, hasting incraasas stability and tands 
to ratard vartical axchanga. Thus axcapt ^ara tha coolad watar ramains at 
tha surfaca and flows aquatorward. soma othar Mchanism. indapandant of tha 
surfaca hasting procass (and in fact countaractlng tha lattar's stabilising 
influanca) must ba prasaiit to vxoosa uaapsr .oyars co tha atmosphara. Only 




s 


th«n it it pottitolt to futl tht wtnn wtttr. attt fonation proeoat at tiM 
xoquirod rttt and allow tha global eyela akatdtad in Fig. 1 to ba ooaplatad. 

Both ataotpharic and ocaanic haat budgatt show that tha priMipal aarhanl an 
aupplying coolar watar to tha turfaea layar ia aquatorial apwalliiig (Haacanrath, IMOf 
1982> Wyrtki, 1981). Of laaaar inportanea it eoaatal npMalling. Cold watar 
advaction in tha aurfaca layar from hi^ to low latitudaa, not involving daaoant 
aaaoeiatad with cooling in tha firat plaea. ia alto aignificant* but eontributaa 
m»dh lata on a global teala than aquatorial upwalling. 

Thua tha main faaturai of tha global ocaanic haat tranaport eyela ara 
known. It ia not claar at all, howavar, what tha eontrola of thia eyela ara. 

What datarminaa tha maaa tranafar rata M , tha tamparaturaa Tp and , 
and tharafora tha nat haat tranaport Q„ ■ M(T^ - T ) 7 Why doaa tha ocaan 
not carry an avan graatar fraction of tha total atmoqpharic-ocaanic polaward 
haat transport than it doaa, whan solar radiation is absorbad in tha first 
instanca in tha top layar of tha ocaan, and whan tha haat capacity of tha ocaan 
is so much graatar than that of tha atmosphara? Tha proximata causa of tha 
ralativaly larga atmospharic haat transport is evaporation from axtansiva 
araas of tha tropical and subtropical ocaan, and th* subsaquant ralaasa of 
latant haat to tha atmosphara, a worldwida procass daacribad in datail in a 
ramarkabla assay by Nalkus (1962) . But why should latant haat transfar ba so 
afficiant? Evaporativa cooling could conceivably daprass saa surface tempera- 
turas to tha point where little haat could ba transferred to the atmosphara, 
sensible or latant. 

In most of tha haat budget and haat transport studies carried out so far 
attention was focussed exclusively on either tha atmosphara or tha ocaan. Tha 
above questions regarding polaward haat transport and its partitioning between 


«tBosph«M and ocaan, howavar, cannot ba anawarad If «ia of tha partnara la 
eonaldarad in »om aanaa Inart. Tha faa atudiaa In %dilch both ataoaphara and 
ocaan %rara takan to ba activa tand to ba ao coaplax aa to obacuxa tha iaportant 
phyalcal controla (a. 9 . tha nuaarlcal modal atudiaa of tha Intaracting atao- 
apharic and ocaanlc boundary layara by Pandolfo and hla aaaociataa, Pandolfo 
and Jacobs, 1972; Brown at al. 1981). Undoubtadly, a fully raaliatlc aiaulation 
of air-sea interaction requires this degree of complexity. However, it is 
likely that considerable insist can be gained from a much sioqpler analytical 
model in tdiiCh sosie of the sore ecaplmx processes are suitably parameterized. 

The present study is an attempt to formulate such an analytical model of the 
Interacting tropical and subtropical atmospheric and oceanic boundary layers, 
explicitly taking into account advaction and tujdaulent transfer processes only, 
and restricting consideration to shallow surface mixed layers. 

In order to avoid explicit consideration of such complex thermodynamic 
processes as cloud fomation, gas and cloud radiation, it is necessary to 
restrict the atmospheric part of the system considered to the mixed layer below 
cloud base. The potential tanperature 6 ^ and specific humidity "above" 

cloud base then become external parameters representing the end product of 
complex atmospheric interactions. This means that some key physical processes 
are at best crudely parameterized by 6 ^^ and q^ , limiting the insight one 
can gain from the model. Nevertheless, the dovfcle mixed layer model yields 
more insight than a single layer one, and supplies some limited ans%rers to 
the questions raised above. 

Although mixed layer models have been widely discussed, a number of pit- 
falls must be avoided in their formulation. Because the literature contains 
many misleading or downright erroneous statements on mixed layer balances, the 


7 


OF POOR QUALITY 


problem is eerofully considered in e generel wey before writing down the 
belenee equations required in the argument. 


Conservation laws with open boundaries 

Consider the atmospheric and oceanic surface sdjced layers in contact* 
bounded above and below by a cloud-base inversion layer ud a diffusiem floor* 
at height z - Z(x,y,t) above sea level and depth z - -h(x*y*t) below 
(Fig. 2) . The equation of continuity and thr conservation law for a 
conservative property in either layer will bj %rritten* using suffix notation 
on this one occasion: 



at 





X) 



( 1 ) 


where x is the property in question and F^ are turbulent flux components. 

The Inversion layer and the mixed layer floor are not fixed either in space 
or relative to the fluid so that Z and h are not* in general* constant, 
and there is* in general, a non-zero velocity of advance relative to the 
air or water: 


H ^ 0 
- II + w(h) 0 


( 2 ) 


>diere w(Z) , w(h) are vertical fluid velocities at the two interfaces 



It is eustoMry to rofor to snd ss "ontrsliuMnt* ^roloeltios la 
saslogy with plvms snd jot ontroinaont problOM* Appliod to s doasity iatorfaeo* 
or 0 soporotion surfoeo dofinod ia oono othor way, tho tom is la iMay rospoets 
o sdsnomr, boeouso it inipllos tho kind of OM<^y progross Msoelotod wltii 
tho growth of jots ond plvMs. It is a»ro illvadnotlng to think of soch oa 
into'rfoco os o shock* front propogotlng rolotiiro to tho fluid in oithor dixoctioa* 
i.o. with positive or nogotivo w^ or w^ . Mewovor* bowing to eustooii tho 
usvAl terminology will bo rotoinod hero* end nogotivo voluos of w^ ond w^ 
will bo lobollod ’'dotroinment* . 

Dopth*intogrotod bolonco oquotions ore used in tho discussion below. 

Given that tho upper ond lower boundaries of the systems connidorod are "open", 
i.o. Z and h vorioblo, care is required in tho integration of tho 
conservation laws, observing relationships of tho kind,* 
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Writing 


0 - f u 
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V ■ I V da 
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( 3 ) 


( 4 ) 


ont rtftdily find* 
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wh*r« is ss dsfinsd in Bq. (2). Ths diqpth intsqrstsd oonssrvstion Isw 
bscomssi 

H * ^ 'f* * ^ 'f" *. X'*’ - ® 

whsrs t is ths div^rqsncs of "diffusivs" horisontsl flux (including tuxbulsnt 
and "shssr** diffusion) i 


z 



o 



ds 


-/ 


V c 
(V - Ji (X - ^ 


ds 


Ths svolution of ths boundary layers will bs suppossd sufficiently slow 
so that t stay be ignored. Alsoi ths layers will be supposed well enough 
mixed to write with a good approximation: 


X 
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u ■ — 


V ■ — 


(S) 


This aj^liss below ths inversion layer. Above that layer the concen* 
tration is different, X * * **y* Given this discontinuity# x(2) in 

Eq. (6) is indeterminate unless one soecifies the sign of w^ . For w^ positive# 
X(2) • Xy » otherwise • X « 
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x(«) ■ x(*)^ ■ 
x(8) - x(*). ■ X 


> 0 ) 
(v^ < 0) 


( 0 ) 


For th« casa of potitiva , Eq. (6) nw xaOueat to 
* « ■ 


( 10 ) 


t^ara 


dt 


lx 

dt 


+ u 


lx 

dx 


♦ V 




and F (Z) and F (0) ara vartical turbulant fluxas across tha Invarsion 
s z 

and tha saa surfaca raspactivaly. A similar traatmant of tha watar sida 
yialdst 


31 + . V 

dx dy w dt 


If ■ 


+ F^(-h) 


( 11 ) 


t^«ra is tha concantratlon balow tha mixad layar floor, 

aqain for tha case of > 0 . A source tarn is 

readily added to tha right hand sida, a.g. to represent enthalpy gain 
by net radiation. Xn tha case of datrainnant. i.a. a convergent mixad layar 
tha aquations remain correct with or sat equal to aero. Physically, 
tha average mixad layar tamparatura is unaffected by a contraction of tha 
system's boundary. It is easily possible to miss this ooint ^Behringer and 
Stonssal, 1981) . 
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Hlx«d l«y«r aod«ls tetcribad in tha litaratura ara oftaa fomu- 

latad by postulating F (Z) to ba aaro (a.g. Tannakas and 

z 

Drladonkt, 1981; Nillar and Kraus « 1977). This is an iqpproxination inappro- 
priate in tha prasant contaxt. At cloud basa, in tha praaanca of paiMtrativa 

convection, tha turbulent flux F (Z) often greatly exceeds tha entrainsiant 

z 

rate “ X) (a.g. Batts, 1976). 

Furthermore, in a nunbar of publications tha "entrainsient flux" 

• X) has bean identified with tha turbiilant flux F^(Z) , a vary 

confusing step after first postulating vanishing F (Z) . "Entrainment flux" 

z 

arises from the choice of a system boundary moving relative to tha fluid; as 
more extraneoiM fluid is incorporated within tha system, tha letter's average 
temperature, or humidity, or whatever scalar property, changes. Rapid 
equalization of tha property within the system implies considerable turbulent 
flux divergences. However, this is perfectly consistent with zero flux at 
the chosen moving boundary: peak downward heat flux in the atmospheric mixed 
layer (to take a concrete example) occurs some distance below the zero flux 
level (see e.g. Ball, 1960). Confusing "entrainment flux" with turbulent 
flvix makes it difficult to reconcile mixed layer or "slab" models with more 
realistic continuum models, a very undesirable outcome (Deardorff and ^iahrt, 
1982) . 


Turbulent fluxes 

Atmospheric fluxes oi interest are those of heat and %«ater vapor. 
Across the air-sea interface these are described by the well-knotm bulk 
relationships (e.g. Brutsaert, 1982); 


Y 


u 




H-(0) 




(12) 



(0) • c 
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«*(q, - q) 


wh«r* R is ssnsibls hsst flux In W * F tits ssm in "kinsMtie* 

S HS 

units, u^ is friction velocity, sstursticn specific husddity st the 
sss surfscs tsnpsrsturs, ^ Stsnton nunbsr, and c^ Dalton nusbsr. Over 
a snail range of tenperature the saturation specific huniditytenperature 
relationship nay be linearixedt 


where 6 is sea surface tesperaturs excess over sons specific reforenoe 
tenperature. F.>r 25*C as reference, q^^ - 0.02 (kg vapor per kg noist air) 
and Y * 1*20 x lo"^ . The atnospheric nixed layer tenperature 6^ is 
to be understood as potential tenperature. 

The turbulent flux of heat and vapor across the inversion layer at 
cloud base can be sinilarly paraneterlsed (Betts, 1976) i 


H. (Z) 

T= w. (0 - 6 ) 

pc * a u 


a pa 


- >'.<■! - V 


(14) 


where 6^ , are potential tenperature and specific hunidity above the 
Inversion and w^ is a nass transfer velocity associated with the "venting" 
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of tho boundary layor by panatratlva eonvoetlon. This proeaaa brinpa about 
an axdwnga of air aczoaa tba intarfaea in addition to any antrainaant* i.ai 
it tranafara haai and aoittura avan whan thara ia no antraiaMnt. 
dz/dt - w(Z) * C . Aa pointad out aarliar# in aixad layar aedala this flux 
cooponant is of tan ignorad, a stop cartainly inadidsaibla at cloud basa. 

Zn tha haat balanca of tha nixad layar in watar tha haat flux 
Ng(0) appaars at on tha air aids, and is in kinaaatic units appropriata to 
tha watar sidat 


® - c u (6 - e ) 

H * a' 0 c 


0 c 

w w 


( 15 ) 


with 0 dansity and tpacific haat of saawatar. 

Tha latant haat lost assoclatad with avaporation isi 


Pw «=w 





( 16 ) 


whara L is latant haat of avaporation. Tha nat haat gain of tha saa 
surfaca by radiation will b« datignatad [k . 

Tuxbulant haat flux across tha ocaanic nixad layar floor is gar.arally 
thought nagligibla, but for contittancy it will ba ratainad for now: 

Fj(-h) • 

whara w^^ is a nixad layar floor nass or haat transfar valocity v-<i 
6^ it tamparatura balow that floor. Tha transfar valocity w^^ nay bacoM 
significant whan tharr is turbulanca on both sides of tha intarfaea, at often 
happens in shallow seas. 
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>%l«nc» •qufctioni OMOMN. mm » 

orNMQUAurv 

Thft balftne* •qoatiom forsnthAlpy in water and in air, and for water 

vapor in air now beeoaw, putting 6 or q for x ^ (3.0) , and eidsetituting 

the various flureei 


a6 K 

(®K • ®> + 7rh~ 

at w *• h pc 

w w 


- - e^) 5*-|i - L 


w w 


w w 


(18) 


d6 


dT ■ <". * '*« - V * 'll “.<* * V 


* * ■ <1u ■ * 'o “‘‘S * ’’ 


It is convenient to introduce the following normalized variables: 


v^ “ c_ u^ 


P. c 

a oa 


* H * j c 
w w 


(velocity scale in water) 


W 4* W 

w ** 


(mass transfer constant in water) 


w + w^ 


B - 




C 

-a- 1 


(mass transfer constant in air) 


(ratio of Dalton and Stanton numbers) 


( 19 ) 
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t - 


s* e_ Y 




OmOINAL MM IS 
or POOS QUALrrv 

(nondlMMlonul h««t) 


q - q. 


(»ix«4 l«y«r v«t bulb tuap^rutur* at tat luval) 


'du 


q^ ~ q^o (wat bulb taaparatura abova cloud baaa* at 
" laval praatura) 


r 0 c V 
w • 


(tantparatura scala of radiant boating) 


Xt indicatad abova. tha ratio V will ba •um>oaad unity and ia not 
ahown in tha aquations balow. Tha raplacamant of aoacifie hunidity by wot 

bulb tamparatura impliaa tha choica of a rafaranca taaiparatura « at which tha 

1 

aaturation apacific humidity ia q (Eq.13). For aimultanaoua validity of 

SO 

all thraa Eqa. (18) it ia. of couraa* nacaaaary to chooaa tha aama rafaranca 
tai^aratura and eonaidar 6 . ate. daparturaa from tha rafaranca atata. 

In contidarinq ataady atata balanca or abort tarm. local changaa it ia 
convaniant to auppoaa tha air tamparatura abova cloud baaa conatant and ^ooaa 
it for tha rafaranca twiparatura. Tha aquationa aimplify aomawhat by writing 
in auch eaaas: 

0 ■ 0 
u 

( 20 ) 

q 5 q (6 ) 

^ao a u 

Tha mora ganaral formulation, with 6^ ratainad aa a forcing variable, 
la raquirad in conaidaring aaaaonal and latitudinal changaa. 
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le 


In tnms of tho noraallsod vnriablos Bqt.(18) bo com t 


h ^ 

dt 


a(6j^ - 0) 


e ♦ 0^ - 1(0 - 0^) + 0^ 


- ae 

Z a 

c u d'fc 

H * 


d0 

Z d 

Cg U. dt 


- V * ® • ®. 

I*(«d„ - V * » - »d 


(U) 


These are inhomogeneous equations for water temperature 0 , potential 

air tenperature 0^ , and wet bulb temperature 0^ . The forcing terms are 

the radiant heating temperature scale 0 , potential temperature 0 and 

%iet bulb temperature 0. above cloud base, as well as water temperature 0. 

Qu r 

below the mixed layer. Kith the coefficients supposed independent of 0 , 

0^ and 0^ the equations are linear. However, at least the mass transfer 
constant 6 is a function of buoyancy flux. i.e. of the air^sea temperat\ire 
difference 0-^0 and of the humidity 0. . In the following, the problem 

A u 

is discussed first for fixed 6 . Later, the dependence of 6 on buoyancy 
flux is considered. Other coefficients in (21) will be supposed constant. 

Steady state solution 

One would expect a time^independent solution of Eqs. (21) to 
model annual average conditions. Such a solution is readily found after 
setting the left-hand side of the three equations zero: 
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mac is 

wot QUALITY 

(1 + 8) (00 «*> e ) *»• 00 -*• 0ft 0 
h r u an 

0(1 •»> 0 ) ♦ 0(1 + ft) 


0 

a 


0 4 00^ 
1+0 


( 22 ) 


0 + 00 . 
du 

1 + 0 


Kay paranatars in thaaa aaprassiona aza tha two notidiiaapsional naas 
transfer constants o and 0 . A significant result follows at once froa 
the definition of o and typical characteristics of subtropical and tjrcqsical 
nixed layers, listed here in Table 1. The scale velocity in water, , is 
of order 0.3 x 10 ^ m s ^ . At the mixed layer floor, tite aiass transfer 
velocity , tdiich parameterizes turbulent flux, is cexrtainly much sneiller 

than v^ . The entrainment velocity w^ is comparable to v^ only in regicms 
of intense upwelling, primarily near the equator. Wyirtki (1981) estimates 
to be 10 ^ m s ^ at the equator, corresponding to a » 3 . Similar upward 
velocities are found in narrow coastal upwelling regions (Hooers et al. 1976) . 
However, over the subtropical anticyclonic oceanic gyres the surface layer is 
convergent, so that there is detrainment and w^ « 0 applies in Eqs.(18), at 
least on an annual average. Thus over most of the tropical and subtropical 
ocean a « 1 • In these locations terms multiplied by a may be dropped 
from Eqs. (22) . 

With the choice of reference state above cloud base (Eq.20), and a • 0 , 
the followinr expression is obtained for sea surface temperature: 


(1 4 jpe 

0 


r 



4 14 



(23) 
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It 


now contains only thzc9 nondiMnsionsl psrsastsrs. 0ns of than* 
it. /8. . is s ratio of two forcing toms, raprosanting raspaetivaly tha 
tasqparatura daprassing affact of avaporation (nota that with tha ehoioa of 
rafaranca stata according to Eq. (20) 6^ is nagativa) and radiant haating. 

As shown in Tabla 1» typical values ara it. « - 17 K and 6 • IS K , 

CIU X 

giving a ratio som«ihat graatar than ona, with a nagativa sign. Tha saa 

surface tanparatura is generally within 1 or 2 K of tha tenparatura above 

cloud base. This iaplias that B cannot be snail cooparad to mity* i.a. 

that w^ + w^ has to be of order c^ u^ . Physically* tha heat transfer 

across tha cloud bast inversion has to be about as efficient as across tha saa 

surface* if naasurad by tha raspactiva haat transfar valocitias. 

Tha affactiva haat transfar velocity w 4- was earlier seen to 

consist of tha inversion layer velocity dZ/dt * the subsidence velocity -w * 

and tha velocity w^ paraaatarising txartnilant axdianga. Of these* tha first 

two can be raadily estimated from data susnarisad by Malkus (1962): dZ/dt 

•4 -1 -4 -1 

no more than about 6 x 10 ms * -w about 3 x 10 ms * or a total 

antrainnant velocity w^ of about 10*“^ m s‘^ . Tha saa surface haat transfer 

•2 *1 

velocity c u is* on tha other hand* typically 10 ms . Xt follotis 

H * 

that tha mass transfar velocity paraaatarising turbulent flux at cloud base* w^ * 
has to be an order of magnitude graatar than w^ . 

Over tha Vanaxuela rain forest, a careful analysis of observations by 
Batts (1976) yielded a mass transfar velocity w^ of 0.13 m s*^ . For c^^ u^ 
as listed in Tabla 1 this would yield B <■ 13 and 6 ■ 6^ • 0 * % 6^^ * 

or mass transfar so affactiva as to erase any property differences across tha 
cloud base layer. However, 'he buoyancy flux in that situation was consider- 
ably larger than found over tie sea surface. The relationship of 6 to the 
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TyplcAl punuMt«ri of trad* wind r«9lon 
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180 

0.01 

0.3 • lo’® 
15 

2.44 • 10® 
1030 

1.2 • 10 “^ 
2.78 
- 6 


20 


9non 0 



buo/aney flux ovar tha ocaan is dlacusaad in graatar Oatail balowi haxa tha 
tantativa conclusion is that 6 is of ordar unity. 

Considar naxt ragions of intansa upvalling whara a is of ordar uni^. 
Tha first of Eqs. (22) nay be rewritten ast 



+ B(1 B)’^ ~ - B(1 + B)"^ (1 + 1) 

o + B(1 + B)"^ (1 + 1) 


(24) 


This gives directly the togierature elevation of the nixed layer over 

deeper water. The three terns in the nunerator are all of the sane order, 

the niddle one negative. With the typical paraneters of Table 1 and a ■ 3, 

B*l> 6 . ■>4K one finds 6 • 6 . about 3 K . The value 9 - 8. ■ 3 K 
n n n 

is exactly %diat was ta)cen to be che nixed layer tenperature elevation by 

Wyrtki (1981) in his analysis of the equatorial oceanic heat budget. 

Although the equilibria solution does not sinulate tenperature and 

hunidity relationships in the interacting nixed layers in a fully realistic 

nanner, with B ■ 1 it reflects the partitioning of the radiant heat gain 

qualitatively correctly. With the left-hand side of the first Eq. (21) set 

equal to zero the only tern representing oceanic heat retention is the 

entrainment term, a (6 - 6) . When ‘a is negligible, the radiant heat gain 

h 


is all balanced by sensible and latent heat loss: with significant a , oceanic 
heat retention competes effectively for sone of the heat gain. In the two 
typical cases discussed before the partitioning is 



Subtropical gyre 
B ■ 1 01 B 0 

Equatorial upwelling region 
B B 1 a B 3 

Radiant heat gain, 3^ 

, 

% 

K 

5 

' 15 

"156 

IS 

166 

Latent heat loss, ^(3-6 ) 

13.24 

88.3 

6.82 

45.5 

Sensible heat loss, 9-6^ 

1.76 

11.7 

- 0.55 

- 3.7 

Oceanic heat retention, a(6-d. ) 

h j 

0 

0 

8.72 

58.2 
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1h« major dlffsronc* betwMn th« subtropical gyras and tha aquatorlal 
\qpwalling ragion is that larga reduction of latent haat loss and its 
replacamant by oceanic haat retention in the divergent siixed layer as the 
most important balance term for the radiant heat gain. This remains 
basically true even idien storage and advection are taken into accounts see 
later discussion. 


Adjustment to equilibrium 

Whether or how closely the steady state solutions discussed above are 
approached in a continuously varying system depends on the rate of adjustment 
to equilibrium. This can be determined by finding the solutions of the homo- 
geneous equations (21) , with the forcing terms set equal to zero, 

6. ■ 6 * 6. * 0 , and still with the reference ten^rature choice of Eq.(20) 

h r du 


0(1 a + i + e’^D) - 0 - H0^ ■ 0 

a d 

0^(1 + B+D)-0«O 


(25) 


0^(1 + B + D)-0»O 
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vh«r« 


D • 



_d_ 

dt 



For the typical values quoted in Ta^^le 1 ■ 278 , so that e is 

small con^ared to unity. Physically, e represents the ratio of the heat 
capacity of a unit area column of air (hei^t Z) to that of water, depth h . 
The detenninant of (25) must vanish for homogeneous solutions to exists 


(1 + B + D) , (1 + B + D) (1 + o + 1 + e“^D) - I 


- 1 - 
« 


(26) 


One root is clearly 


« - (1 + B) 


(27) 


Further inspection reveals that one of the remaining roots is of order 
unity, the other of order e . To order e these are: 


Dj • - (1 S) - e ^ * g + O(e^) 
Dj • - e(o + B £" “-| ) + O(e^) 


( 28 ) 
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Su^ws* ROW thot th« oyttwi ttart* froR mm «rbttr«cy non-oqaili)»riM 

•tat* at tlM t ■ 0 » glvan by initial taaparaturaa 6 , and 6. , 

O mO QO 

Taapantura and husddlty changaa in tha eooraa of adjuatmnt to aquillbriuR 
ara than givan by 


8*^ - " ^ji •'i^^ 






( 2 ») 


• H. 



whara an ovarbar dacignatas tha praviously found staady stata solution, 
j is succassivaly no subscript (watar) , subscript a or d , and i ■ 1,2,3, 
with sumnation ovar 1 impliad. Aftar soma routlna calculations ona finds 
tha tamparatura parturbations 6'^ axprassad in tarms of thair initial 
valuas , to saroth ordar in e < 


.. . '• * .-V, t . 

« , •‘-TTT'* ' * rrr* ttttt 


1 -’ 


kjt 


1 ■* B • 


- c V * ' C. C . 


C 1 -k t 


d 1 + i 


1 + i 1 + A 1 + Bi 


(30) 


C 




6' - C a 


-)C3t 


£J*Q**'*ai mkk m 
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Vh* SM turfaet twiptntur* ■tiiply adjust! to Its s^UteltM iralus on 

ths slow tias sesls . Ths sir tsapsrsturs and spaeifle hwldity initially 

also chanQs on tha fast tint sealas and k^”^ » but at t » k^*^ thsy 

just follow tha saa tanparatura. For 0*0, 6*1 and data in fabla 1 

tha valua of k^*^ is 0.81 x 10^ s » or 102 days. Tha tisa soalas k^*^ 

_2 

and k^ ara tha saaa to saroth ordar in c and do net dapand on a For 
6 ■ 1 and data of Tabla 1 thair typical valua is ■ 3 x 10^ s or 8,3 hrs. 

Physically! diangaa on tins seals k^~^ raprasant tha adjustnant of tha 

•X 

spaeifle h\anidity in tha atmospharic stlxad layar, thosa on seals k^ of tha 
air tas^xmtura, «d\ila tha procass on tins seals is elaarly tha adjustnant 

of tha hast contant of the ocaanlc nixad layer. 

Tha above calculations have ravaalad tha tlM sealas of adjustnant* but 
ara realistic models only at a coastline (a. 9 . tha kfrican coast of tha North 
Atlantic) where tha air blowing from land adjusts to tiia balances dictated by 
tha oceanic mixed layer. Xn mid-ocaanie regions local changes 0/2t) arise 
in response to diurnal and seasonal variations of radiant heating* as wall as 
of the atmospharic tanparatura and huaddity above cloud base. Thasa nay be 
modeled by harmonic forcing tarns, Advactive ^langas result from the global 
scale gradients of tha forcing terns, coupled with along*gradiant flow. The 
tine dependent response to general* variable forcing is also readily written 
down and serves as a basis for considering local effects of cold air and water 


advaction. Thasa effects are discussed in tha next two sections 
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Rtiponit to PTlodic forcing 


^ITY 


-1 -1 

Th« typical tiJM scales and ara fortuitously close to 

the two principal frequencies of radiant heating » the diurnal and the annual t 


- Uy ■ “ 1“’* 


0) - — ■ 2 X lo"’ s“^ 
A year 


Tc consider the response to periodic forcing at these frequencies let 
the radiant heat gain be written 


e 


r 


6 e ' iwt 

r re 


(31) 


It is only necessary to calculate the response to tht periodic part, setting 

at first 6. ■ 6. ■ 6 *0 in Eq.(21), retaining for now the reference 

du n r 

temperature choice of Eq.(20) and supposing 6 , 6^ and 6^ to vary 

periodically: 


e 


iu)t 


e 

a 




(32) 


where the amplitudes 6* , etc. may be complex. Equations (21) reduce tot 


(1 + 6 + io) - 6' 

6' (1 + 6 + io) ■ 9’ 

d 

e* (1 + a + 1 + e“^ io) • e* + le* + e* 


( 33 ) 




0 ■ 


u z 
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(34) 


is ths nondlncntional valus of th« forcing iroqucncy. Th« solution is 




6' - (1 + S + ia) 0' 


(35) 


v)h«rc m , ^ ars raal, dafinad by 


n a^^ ■ 0(1 + i) + ad + 8) - + 


-1 


(36) 


+ io[c (1 + 6) 1 o ♦ 11 


Suppose that tha nondimansional fraquaney o is of ordar unity: 
this is tha casa for diurnal forcing. Tha loading tarns in Eq. (36) ara than 


n a^*** S - e"^ 0^ + e‘^(l + 8) 


io 


(37) 


so that 


-tan 


"IH^I 


-1 2 
BBC O 


1 + 


(1 + 8 ) 


2l 


.2 j 


?or tha typical valuas of Table 1« is 4.63. With 8 ■ 1 > the 

corresponding phase angle ^ is -23.4*, and according to Eq.OS), 0* and 3* 

A a 
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or mu 


sm 

WWLfTf 


l«ftd 6'^ . howttv«r, th« r««ponM !• of nogllgihl* agplitudo# proportional 

-1 -4 

to m , which !• typically of ordar 10 

In tha casa of saaaonal haating (o of ordar e) it ia not raalistic 

to Bupposa 6^ - constant t air teaiparatura also variaa significantly abova 

the inversion, in phase with 6^ . Thus although ■ 0 nay still be 

chosen for a reference tanperature, 6'^ must be supposed variable, and so 

must 6' . A more complete version of Eqs. (33) is, neglecting now 

GlU 

quantities of order O t 


e* (1 + B) ■ e* + se* 

a u 


e; (1 + 6) - 6’ + 661 

d du 


38) 


6* (1 + a + 1 + e~^ io) 


e* - £0’ - 6 ' 

a d r 


Under this set of idealizations the sea surface tenpwrature is found 


to 


S' ' i's ♦ s> % * ewiJ 

IQ € L ^ 


(3V) 


with 


ine^^ *6(1 + £)+ a(l + 6)+io e”^(l + 6) 


or ^ » tan 


~li_e^(l + 6) 


^6(1 + 1) + 0(1 + 6) 


a^d + 6)^ + 16(1 + £) + 0(1 + 6)]^ 


Quality 


For «nnual torcinv ^ * 0.012 , «nd tht «bov« r«lation»hipt givt 

for tho physical paranattra in Table 1« with a ■ 0 » 0 ■ 1 i 


0 - 60* 
n • 7.67 

.2 

Tha typical amplltudia of tha annual forcing la 75 N n , corraaponding 

to 6* ■ 6.25 K , whlla 6* la of order 2 X , 6' tha aaaa. Tha 
r u QU 

raaponaa aatplituda 6' la than for 0 ■ 1 2.7 X . lagging behind tha forcing 
by 60* or two months. Tha first two of Eqs.(38) yield 6* ■ 0*^ ■ 2.35 X, both 
in phase with tha sea surface temperature. 

A quantity of soma practical interest is tha seasonal heat storage in 
tha oceanic mixed layer, represented in Eq.(21) by the term 


s 


h de 


<40) 


In the above periodic solution this varies as i<76' , i.e. it leads 
the forcing by (V2 - 0) .or typically 30* , corraspondinq to a month. 
Purthazmoro . at this (annual) frequency the storage term is only of order 
unity in the first of Eqs. (21), i.e. in the heat balance of tho oceanic mixed 
layer, not in the heat or vapor balance of the atmospheric layer. With the 
typical amplitudes used above, tho storage term amplitude a 6* 
corresponds to a dimensional value of about 90 W m ^ . The seasonal cycle 
of temperature changes and heat storage is clearly of significant amplitude. 
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Cold waf r tad >lr advction 

Advoctivo changt* vay b« troAtod mudi u local onaa« by auppoaing that 
tha forcing tama 6^ , 6^^ and 6^ ara functions of location* and raplacing 

tina by spaca darivationsi 


__Z d . 2V d . d 

Ch dt ■ Cjj u. dy = dY 


K D 


h _d_ 
dt 


hv _d_ . d 

dy * <K 


< 


hv Cjj Uj, 
2 V 


2 


(41) 


Hara tha gradiants hava for sin^licity baan supposad to point along 
tha y axis* vdiich is tha diraction of positiva advaction valocitias V and v * 
in air and watar raspactivaly. Tha forcing tarns in Eqs.(21) 6^ , 9^ and 

6^^ will all ba supposad functions of tha coordinata y alona: 

0(l+i+<D)-0 - 10. • 0 (y) 

A or 


0^(1 + $ + D) - 0 - S e (y) (42) 

A U 

0^(1 + $ + D) - 0 - S0j^j(y) 

(o S 0) 

The case of negligible upwalling (a ■ 0) has baan assunad for 
simplicity, because later discussion of advaction will deal mainly with tha 
subtropical gyres. Tha typical value of tha parameter K = (v/V)e ^ is now 
not large, but of order unity, on account of tha small typical value of tha 
advaction ratio v/V . Usual advaction velocities ara v ■ 0.05 m s ^ and 


K 


i 
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V « 5 n , giving ic ■ 2.78 . Th« nondiMnsional dlstane* Y is seslsd 
by ZV/e , which is typically 300 ka. 

Upon elininating 6^ and 6^ from Eqs. (42) ons finds tha single 
aquation for saa surfaca taxnparatura: 

ic e + [1 + I + <(i + B)i D e + B(i + A) e - ♦(y) (43) 


tdiara ^(Y) is the affective forcing function for saa surface tasparature: 


(^{Y) - (1 + B D) e «*• ee + b a e 

r u ou 


The characteristic equation derived from (43) has roots r^ , r^ 
given by 


r, + r - - ' ^ — +1+6’ 

12 ic 


(44) 


^2 


(1 + A) 


The solution of (43) is: 


6 


•c(r2-rj^) 




4>(Y - n) dn 


(45) 


0 

The result shows that the forcing terms over a "badeward** (Y* < Y) 
sector of the y axis, from where the advection comes, influence tha local 
temperature. With r^ , r^ of order unity, the width of the influence zona 
is the scale of the Y vari^U^le, earlier seen to be typically 300 )cm. Ttds 
is small on a global scale and it is therefore realistic to replace the 
forcing function 4>(Y) by its linear expansion: 
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♦(V - n) : ♦(») - 1 ^ 




(46) 


Substitution into (45) then yields t 


e - 




»(Y) 

1C r,r_ K 2 2 dY 

12 


(1 + 6) 6 + ee + B 1 6. , 

r u qu ^ 1 


d6 


S(1 + 1 ) 


B(1 + 1) dY 


(47) 



1 S 1 d6 
S(1 + 1) ' dY 


The first term on the right is exactly the equilibrium solution 6 
(Eq.22) for a * 0 . The second and third terms represent a temperature 
perturbation associated with cold water and air advection, \dtich arises idien 
the forcing functions vary in the direction of advection. 

Further calculations using Eqs. (42) yield for the dry and wet bulb 
temperatures : 


e 

a 


+ £ 9 ^ 

U 


1 





0 


-r n)^ 

e U»(Y - n) dn 


where 4i(Y) ■ (1 + £) 6 + 6(1 + £ + <D) (6 + £ 6 , ) 

r u qu 


r I -d + B) 

e - e . - B ' 8 (Y - n) - 0 . (y - n) i e ari 

a u , u qu 


(48 
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or, using a llnaar sxpantlon again, on tha sedol of Bq.(46)t 


6 + Jl 6 - 

a d 


if f 1 , , . 1 4 B 1 d» 
B(1 + L) “[6 B(1 4 L) J dJf 


(49) 


6 

a 


®d- 


1 4 B 


(6 


®du^ 


B 


(1 4 B)' 



d 6 


du 


d y 


It is raadily verified that the last expressions also consist of the 
equilibrium solution plus terns proportional to the gradients of the forcing 
functions. In case the latter are not known very well, the perturbations can 
also be expressed directly from Eqs.(42), in terms of the gradients of the 
dependent variables themselves: 


e - e-e* 


6 

a 


6 

a 


D 6 4 0' 

a 

~rrs 


(50) 


6 


d 



06.46' 

U 

14 8 


%«here 6' is the sea surface tenperat\ire depression: 


6* 


ica 4 B) 0 04 D & 4 I D 8 . 

a g 

B(i 4 i) 


While the gradients of the forcing terms are hard to estimate 
there is good evidence on the variation of sea surface temperature in tlie 
direction of the likely surface advection over at least part of the North 
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Atlantic •ubtroplcal gyra. Ovar tha northarn half of this gyra surfaea drift 
is nora or lass along and "up" tha surfaca tanparatura gradiant. According to 
tha diarts of Bohnacka (1938) that gradiant has a typical nagnltuda of 
2K/1000 ka (2 X 10 K a ) , so that using aarllar astiaatas 


Tha air-saa tanparatura diffaranca and tha dry bulb— wat bulb diffaranca 
do not change vary much, and one may supposa in a first approximation t 


d 6 

a 

dY 


de. 


dY 


u 

dY 


(51) 


Tha typical magnitudes of tha advaction tarms in Eqs. (42) than become 

K D 6 ■ 1.67 K 

D0 ■D0. "O.SK 
a d 

With the aid of Eqs. (50) these give the following corrections to tha 
equilibrium temperature: 

0' - - 1.51 K 

0 ' - 0.'- - 1.05 K 

a a 


showing that cold water and cold air advaction tends to reduce the sea-air 
tenperature difference , because it causes a greater tenperature depression 
in the water than in the air. 
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Radiation h««t lo«« of ttv atiaotphTlc nixmi l«y#r 


Ona obsarvatlon i« difficult to raconcila with tha haat 

balance equation for the atmospheric mixed layer as formulated so far 

(second of Eqs.21 or 42). This is the fact that 6^ is almost always less 

than either 6 or 6 . Given efficient heat transfer across both interfaces, 

u 

the right hand side of the second Eq.(21) is typically greater by 1 or 2 K 
than can be explained by cold air advection. One is led to consider heat loss 
by gas radiation H , another potentially important forcing term in Eq.(18). 

AIT 

In a normalized form appropriate for inclusion in tha second of Eq. (21) this 
term is: 


e 

ar 


ar 


0 c c u 
•^a pa H • 


(52) 


The magnitude of 6^^ may be estimated from data sumnarized by Fleagle 
and Businger (1963) , according to which the air at low levels cools ]t>y radiation 
typical y at the rate of 3 K per day or d9/dt «3xlO^Ks^ . The corre- 
sponding heat loss term in Eq. (18) would be 


ar 

P c 

a pa 


- Z II - 1.8 


X 10 ^ K m 


-1 


With c^ u^ ■ 0.01 m s ^ this gives a typical value for the normalized 
variable 6 of - 1.8 K or more or less the magnitude required to balance 

AIT 

the second Eg. (21) . 

The modification of the equilibrium solution (Eq.22) on account of the 

6 term consists of the inclusion of 9 in two of the three numerators, 
ar ar 

giving correction terms of 
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6 ' 



8(1 + A) 


e • 

a 


e 


1 + 8 


S?'22£ » 
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(53) 


Typical valuta art - 0.9 K for 6^' , - 0.5 K for 6* , 

con^rable to the effect of cold water and air advection^ and contributing 
further to the sea surface tenperature depression, vihile also modifying the 
sea-air temperature difference, this time in a positive direction. 

Heat and mass transfer through penetrative convection 

In order to complete the analytical argument, it Is necessary to consider 
the physical processes responsible for determining the value of the mass 
transfer coefficient 8 , and show that the order one typical values supposed 
above are plausible. The problem of mixed layer deepening (in the atmospheric 
or oceanic application) has an extensive literature (e.g. Niiler and Kraus, 
1977; Tenne)ces and Driedonks, 1981) but it is almost exclusively aimed at 
density interfaces across which the turbulent flux is negligible, and entrain- 
ment is the only mechanism of heat and mass transfer. As pointed out above, 
this is not the case at cloud base, and a different approach must be adopted. 

Given the intense penetrative convection at cloud base, it is 
reasonable to suppose that the vertical motions that "vent" the mixed layer 
have the velocity scale of free convection (Deardorff , 1974) : 

w^ - (BZ) (54) 

c 


vdiere B is the sea surface buoyancy flux, due to the transfer of sensible 


h«at and also of watsr vapor (Brutsasrt, 1982) t 
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B 


c u 

H * 



(55) 


with T tha rafaranca abaoluta taniparatura. 

nia slnplest paranatarlzution scbama for tha nat mass transfar 
coafflciant at cloud basa is: 


w + w . ■ X w 
a • c 


(56) 


whera X ■ constant. This cartainly doas not do full justica to tha 
physics of panatrative convection » tha nat mass transfar baing in addition 
dapandant on the distance Z from tha lower boundary, tha buoyancy saltus 
b at tha top of tha mixad layer, and possibly some* other factors, i.a.: 

2 ^ 
w 

X ■ func (57) 

Typically, however, the cloud base inversion is weak and it is not too 

unreasonable to suppose that turbulent transfer processes in its neighborhood 

aure determined by tha parameters characterizing the free convection regime, 

a supposition that leads directly to Eq.(56) . This amounts to supposing 
2 

that bZ/w is suitably small in tha parameter range of interest, 
c 

In the case analyzed by Betts (1976) in detail the characteristic convection 

velocity was w 1.85 m s'^ , while the nat mass transfer coefficient was 
c 

empirically detemdned tobe w +w. «0.13ms^. This gives X ■ 0.07 , 
a value that will be adopted in the following. 

Combining Eqs. (54) to (56), expressing them in terms of the variables of 
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Bq:(19), *hd substituting ths sguiiibriun rslstionships (22), ons finds ths 
follotring sxprsssion for ths cosffieisnt 0 (rsfsrsnes tsapsrstuxs choies 
secording to Eq. (20) t 


0 ** 6 , 


6^ ♦ 0.61 YT 





^ALny 


IS 


(58) 


whsrs 6^ is s tsinpsrsturs seals t 


e 


* 




T 


(59) 


and 0^ is givsn by (22). 

With ths abovs chossn valus of X ■ 0.07 , T ■ 300 K , and ths typical 
quantitiss of Tabls 1, ons calculatss 

0^ - 0.167 K 

whils 0.61 y T > 0.22 . Ths valus of 6 may now bs :,stsrminsd (supposing 
stsady stats conditions) from So. ($8) and (22). as a function of ths cosffi- 
cisnt o , and ths forcing paramstsrs 0^ » 6^^ and 6^^ . 

With ths typical quantitiss of Tabls 1. and for a ■ 0 ons finds 

0 - 2.3 

Physically, such a rslativsly low valus of 6 (cong>arsd to ths ovsr<>land 

cass of Batts. 1976) is dus to ths fact that most of ths hsat is transfsrrsd 

as latsnt hsat. which causss only littls buoyancy flux. 

For a ■ 4.02 . 6 drops to tsro: at this watsr sntrainmsnt rats. 

with 6 • - 4 K . thsrs is zsro buoyancy flux, and by ths hypothssis of 
h 

£q>(54). vani'ihing sntrainmsnt of air from abovs ths atmosphsric mixsd laysr. 
This asynptotic cass is characterized also by 


(60) 


e 


e. • e. 
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or an atmosf^arlc mixad layar of tha aana tanq^atura aa tha aaa surfaca, 
saturatad with watar vapor, claarly a soaiawhat onraalistlc raault. A 
valua of 6 ■ 1 la raachad at a • 2.57 . lha praclaa nuatarical valuaa 
ara not aignificant in virtua of tha uncartainty of X , which waa daducad 
fr«n a ain 9 la caaa obaarvad ovar land. Howavar, tha calculatlona ahow that 
an ordar ona valua of 6 ia compatihla with raaliatlcally aaaumad propartiaa 
of penatrativa convaction ovar tha ocaan. lha valua of 6 ia kapt nodaat, 
bacauaa thv. buoyancy flux ia low at tha uaual ralativaly high rata of 
avaporatlon. 


Coapariaon with obaarvation 

So far tha anphaaia in tha varioua partial compariaona batwaan analytical 
modal and obaarvation haa baan on tha datamination of tha ordar of magnituda 
of tha diffarant tarma in tha haat and vapor balancaa, Eqa.(21) , tha tranafar 
coafficiants, cold watar advaction, ate. Aftar tha kay paramatara hava baan 
choaan, a mora data! lad axanination of tha obaarvad t««paratura and humidity 
ralationahipa bacomaa poaaibla and ahould anawar quaationa ralating to warm 
watar maaa formation. Datailad ralavant avidanca ia aval labia in tha aubtropical 
trada wind bait, both in tha form of atmoapharic tamparatura and hum'.dity 
profilaa, and tima hiatoriaa of monthly avaraga aaa to air fluxaa, air and 
watar tamparaturaa . 

Tanparatura and humidity profilaa in tha aubtropical marina boundary 
layar hava baan publithad by a numbar of invaatigatora, including Bunkar at al. 
(1949) , takan just north of Puartu Rico in Apral, Pig. 3. Tha followinq 




3 » 


twBp«ratur«s nay ba axtractad from thosa data (rafaranca 6^ ■ 0) t 


e. • - 11.5 K 
du 

e. - - 7.5 K 

Q 

e • - 1.5 K 

- - 2 K 


Tha wind apaad was about 


8 in t 


-1 


from which ona attinatas 


Si 


■ 0.02 n 


-1 

s 


Ttie nat haat gain from radiation, attimatad from climatological data for this 

-2 

time of yaar 200 W m (Bunker, 1976) , rasults in: 

0 • 8 K 
r 


Tha absoluta potantial tamparatura abova cloud basa may ba takan to 
hava bean 301 K . Substituting into Eqs.(21) ona finds: 


& . 9.2 

dt 


_ de 

-5 -T— -0.5 + 28+0 

ar 


d0 

2 d 


u, dt 


- 48 + 3.5 


where radiation heat loss has also bean explicitly included. 
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Allowing also for cold air advaetion as par aarliar astlsiata (d6 ■ 0.3 K, 

for tha high c„ ) , tha sacond and third aquations cons closa to balanea 

with 6 ■ 0.8 . Tha first aquation, howavar, shows that tha watar oolvaon 

was cartainly not in tharmal aquilibriun: supposing tha usual waak rata of 

cold watar advaetion eharactaristic for this location, tha watar column was 

-2 

baing locally coolad at a rata of soma 250 W m , whlla in April a haatinq 

rata of that magnituda is usual. Ona raasen for tha discrapancy is tha 

smallness of ^ 2 ’ ' * diract consaquanca of tha high wind spaad, tha value of 

e u^ baing twice the typical value listed in Table 1. Another reason is 

the large negative value of 9^^ , i.e. the dryness of the air descending 

into the mixed layer. The exanple demonstrates that instantaneous or short 

term average temperature and flux conditions cannot be used to infer tha 

monthly average oceanic heat retention rata. One notes also that the heat 

-2 

gain estimate (200 W m ) on the basis of a monthly mean may not have bean 
appropriate . 

The equilibrium sea surface ten^erature one calculates from (22) , 
with allowance for cold water advaetion and radiation loss from the air, is 
6 > - 4.3 K or some 5 K lower than observed. Given an anomaly of 5 X , 
one expects the normalized cooling rate to be 

- It ■ - h, A0 ■ 1.89 A0 - - 9.5 X 
V* v^ 3 

having used Eqs. (28) and (29). The result agrees almost exactly %d.th 
the first of the balance equations (21) , with observed tc’tg>eratures 
substituted, see eibove. The agreement corroborates the climatologically 
estimated rate of cooling of 200 V7 m ^ . 


OfflOINAL PAQC IS 
OF POOR o’.'^*,nrY 


In spit* of th* thosnsl dis*quilibriom of th* oeoanie mimd l*y*r, 

aumy aspects of th* data in Fig. 3 ar* typical for th* w*st*m North 

Atlantic trad* wind region. Th* s*a surface t*iqp*ratur<i Is here wailly 

slightly higher (0.5 to IK) than the siixed layer air tesperature, but 

lower by about 1-2 K than the air above cloud base: there is usually a 

weak inversion at cloud base. Th* very large value of 6^^ seeiw to 

be connected with th* height of the main trad* inversion. As may be 

seen from th* figure, th* local wet bulb— dry bulb tenperature nearly 

vanishes at the top of th* cloud layer. Correspondingly, the specific 

humidity within th* cloud layer is near th* saturation value at cloud 

top, i.e. at a height in this case of 2300 m. Cloud base is, of course, 

at th* lifting condensation level, or where th* sea level air becomes 

saturated. Th* difference in saturation specific humidities between 

cloud base and cloud top, and therefor* 6^^ , is a function of cloud 

layer depth. As Malkus (1962) discusses in detail, this depth grows 

from the African r.«st westward over th* North Atlantic. The typical 

value of 6^^ ■ - 6 K oi' Table 1 was chosen to represent an intermediate 

stage. Even at Puerto Rico the large wet bulb tenperature depression 

seems excessive: 6. ■ - 9 K should be more typical. 

du 
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Accepting 6. ■ - 9 K for the Puerto Rico locetion* using 6 ■ 1 > but 

GlU 

the more typical air-eea transfer coefficient c u^ ■ 0.01 m s*^ and an annual 
net radiation gain 180 W m ^ , one calculates the following equilibria 
temperatures from Eq. (22) t 

e « 1.32 K 

e - 0.66 K 

a 

6^ - - 3.84 K 
a 

These are clearly wrong: air below clcud base should be colder than above. 
Correcting for radiation loss (Eq.53) one finds the slightly more realistic 
results: 

6 ■ 0.84 K 

6 - - 0.48 K 

a 

9 ■ - 4.08 K 

d 

The calculated sea surface temperature is still much too high. However, 
after allowing for cold water and air advection according to Eq.(Sl) one 
arrives at: 

e • - 0.67 K 

e - - 1.54 K 

a 

6^ - - 5.15 K 

U 

The sea-air temperature difference is now within its typical range of 
0.5 - 1.0 K ; the dry bulb — wet bulb temperature difference is 3.6 K, which 
is somewhat low but not atypical, and the cloud base inversion strength is 
1.54 K or pretty much as observed. Better agreemnt can hardly be expected. 
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Th« aain point it that without allowing for atmotpharic radiation lota and 
cold watar and air advaction tha longotaxm naan tamparatura-hunidity 
ralationahipa cannot ba at all raaliatically ainulatad. 

It it of intaraat to oonaidar alto tha partitioning of tha total haat 
gain 6^ of tha ocaanic mixad layar into aantibla and latant haat lota and 
nat haat ratalnad (tha lattar owing to cold watar advaction) . Attuning that 
tha jutt calculatad valuta realistically rapratant annual avaraga conditiont 
ona finds tha fol lowing: 


(K) % 


Heat gain, 9^ 


15 

100 

Latent heat loss, i(9 

1 

CD 

12.45 

83 

Sensible heat loss 9 

• 

CD 

.87 

6 

Cold water advection 

vh/v^ d9/dy 

1.67 

11 


As remarked earlier, these differ little from what one deduces front the 
equilibrium solution. The Bowen ratio (6 - 9^)/8.{0 -• 6^) is 7%, or a typical 
value. The dominant* of the latent heat loss is seen to be caused by tha high 
value of the material constant I , and the high absolute value of 9^^ (which, 
with 6 of order unity, translates into a large dry bulb— wet bulb difference 
in the mixed layar) . As pointed out above, tha high absolute value of 0^^ 
derives from low saturation specific humidity at the height of the trade inversion 
and is, in a sense, another property of the water substance, although it is also 
an outcome of boundary layer processes determining the thickness of the cloud layer. 

Monthly and yearly average temperatures and fluxes have also been deteirmined 
from sea surface and satellite data by a number of investigators. In the North 
Atlantic trade wind belt a typical location is 1500 km northeast of Puerto Rico , 
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Th« partitioning of th« hast gain lax 


(K) % 

nat haat gain 6^ 17 100 

latant haat loss i(6 - 6^) 14 82 

sansibla haat loss 6-6 0.5 3 

u 

cold watar advaction 6 ■ vh/v. d6/dy 2.5 15 

n ^ 


These data are vary similar to the ones just calculated from the steady 
state model under typical conditions. The only difference of any note is 
a somewhat smaller sensible heat loss, due to a small sea-air teng>arature 
difference. The correspor.dence establishes that the equilibrium solution of 
the boundary layer model, modified by an allowance for cold water and cold 
air advection (and with the inclusion of radiation heat loss in air) realistically 
simulates annual average conditions. 

A perhaps more stringent test of the simple model is a comparison of 

results calculated for annual harmonic forcing with the observed variation of 

monthly mean temperatures. For this comparison it is plausible to take the 

value c^ u^ - 0.008 m s ^ deduced from the annual average temperature- flux 

relationships and use 8 ■ 1 • The amplitude of the forcin<, function is 
•2 

• 75 W m , directly from Fig. 4. Less certain estimates are mixed layer 
depths of Z *> 600 m and h * 50 m , and seasonal fluctuation amplitudes of 
dry and wet bulb temperature above cloud base, 0 ' ■ 6 ' ■ 2 K . The non- 

dimensional annual frequency is with the above value of u^ equal to 
O ■ 0.015 . From Eq. (39) one Tinds: 

4> - 65.62 
m - 9.16 
<>•- 2.46 K 

O e'- 10.27 t; 


PRECEDING PAGE BLANK NOT FILMED 


46 


original * 
OF POOR QUALITY 


$»aA4 


This rsprsssnts s phsst lag of s littls ovor two nonths, a soa surfaco 

t«Bg>«ratura ratponsa anplltuda of 2.5 X and a haat storaga ao^lituda of 

ICJ w . Tha lattar la to ba coagtarad with tha a^plltuda of ocaanic haat 

ratantion, H in Fier. 4. Tha storaga should load tha haating rata by about 
n 

ona month. Tha modal pradictions ara vary closa to tha obsarvad rasults* 
considering tha crudenats c I soma of tha astimatas that want into tha calcu- 
lations. All of the evideaca confirms that 6 •• 1 is a raasonabla astimata 
of tha cloud base mass transfer coefficient. Tha order of magnitude of this 
constant may therefore ba taken as firmly established. 


^^o6es of oceanic haat retention 

What does the above analysis imply about warm water mass formation? 

The process takes place when and where the ocean retains a portion of its 

radiant heat gain, i.e. when there is net heat retention H on top of %diat 

n 

is transferred to the atmospheric mixed layer: 

”n • ”r ■ ^ 'pa‘“ ’ V ' ^ ~ 

According to Eg. (18) the net heat retention is also 

“n ae 

» If ® V <“> 

w w 

The right hand side of this equation shows that heat retention by tha 
oceanic mixed layer may manifest itself in three different ways: (1) as 

local heating, i.e. storage, h 36/3t ; (2) as horizontal cold water advection, 

or haating the fluid as it moves along, h u • V 6 ; (3) as vertical 
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adv«ction, l.«. haating of fluid antrainad from balow, acconodatad in tha 

divargant surface layer. In tha previous section storage and cold water 

advection ware considered in detail in the North Atlantic subtropical gyre. 

Fig. 4 shows net annual heat retention of order 30 W n ^ . According to 

the charts of Bunker and Worthington (1976) or of Hastenrath (1980) , a similar 

net gain diaracterizes large areas of ti.e trade wind Iselt, i.e. the subtropical 

oceanic gyres. However, these values of H were estimated from meteorological 

n 

and satellite data as a difference between positive and negative terns in an 

energy budget like Eg. (18), with f\irther split into incoming and outgoing 

radiation. It is readily shown that the typical error of such calculations is 

-2 

of the same order as the calculated net heat gain, i.e. 30 W m (Weare et al., 

1981) . It is therefore important to examine the oceanic evidence to see vdiether 

the net heat gain values can be substantiated. 

The comparison of the botuidary layer model with observation in the last 

section showed that the temperature- flux relationships near Puerto Rico and 

northwestward can only be understood if cold water advection is taken into 

account, at about the rate corresponding to the meteorological evidence. This 

rate was inferred in the section dealing with cold water advection from observed 

ocean surface tenperature gradients and the generally accepted magnitude and 

direction of the surface drift over the northern portions of the subtropical 

gyre. Other estimates of heat gain associated with cold water advection have 

been made by Stommel (1979), Clarke (1979) and Behringer and Stommel (1981) at 

relatively low latitudes (20*N and 9^N respectively). All these estimates 

-2 

yielded values of order 10 W m , in rough agreement with the charts of 
Bunker and Worthington for the region in question. 

Another way to examine the oceanic heat retention is to calculate the 
global heat balance of the entire mixed layer over a subtropical gyre. 


48 


OmOlNAL PAGE 18 


Putting ■ 0 in Eq. (62) for th« cm* of dotrainnuit and draping 

th* storage tern in a calculation of th* long term mean heat galUf the equation 
is integrated over the area A of a gyre, enclosed by perimeter C. On* finds 
by the divergence theorem and the wat*r>sid* equivalent of Eq.(5)t 



6 dx dy 


(63) 


The second term on the right arises from the horizontal divergence of 
the flow, Eq. (5), an integral of iidiich, adapted to the water side yields: 


r 


9 


c 


h u • 


V 

» 

n ds - ' i w dx dy « 0 

- ) w 

A 


(64) 


In these expressions ds is a line element of the perimeter C , n its 
outward unit normal. The integral balances (63) and (64) can be written down 
^ priori. 

Over the subtropical gyres w is negative and there is net inflow at 

w 

the perimeter. Let that portion of the perimeter with inflow (u • n < 0) be 
designated , the remainder • Inflow and outflow volumes are then 


^ h u • n ds 


Vj ■ ^ 


h u • n ds 


(65) 
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Writln 9 \ w^lght^d «v«rag« t«np*r«t\ir«s of 

nixod layor inflow, outflow and datrainntnt, ona «ay than axpraaa Bq. (63) as 



H 




dx 


dy - [ 





( 66 ) 


This balanca contains bulk quantitias not too difficult to astisiata from 

axisting information. Maar tha southarn boundary of tha svibtropical gyra 

tam^aratura diffarancas ara small. It is not unraasonabla to supposa that any 

wind drift antaring from tha south is datrainad in tha intarior at about tha 

sama tamparatura as it antarad, so that it makas littla contribution to Eq.(66) . 

Total inflow into tha mixad layar at tha northarn parinatar (naar 40*N) dua to 

wind drift a Iona is of ordar 10^ m^ s~^ . niis is haatad, according to 

Bohnacka's (1938) charts, by soma 5 K bafora baing datrainad. Tha corrasponding 

hast gain, distributad ovar an araa of 2000 km x SOOO km (about tha araa of tha 

-2 

gyra) amounts to 20 W m Additional cold watar advaction by southward 
gaostrophic flow incraasas this astimata somawhat, as doas an allowanca for 
cold surfaca watar antaring tha gyra's mixad layar from tha upwalling ragion 
at tha aastarn boundary. Cruda as tha ocaanographic astimata nmy ba, it is 
robust in tha sansa that naithar tha total southward transport across 40*M , nor 
tha tamparatura risa batwaan (say) 40*K and 20*N. nor tha total araa of haat 
gain can ba sariously in arror. Ona can tharafora accapt at laast tha global 
valua of tha nat ocaanic haat gain ovar tha North Atlantic subtropical gyra as 
baing of ordar 0.3 x 10^^ w , as inpliad by tha charts of Bunkar and 


Worthington (1976). 


In th« arM of oquatorial upwolling an important haat rataation machanim 
la tha wanning of watar an taring tha mixad layar from balow, which ia acoome* 
datad in tha divargant aurfaea layar. In Eq.(lB) tha corraaponding haat gain 
waa axpraaaad aa ^ turbulant tranafar valoeity w^ baing 

praanmably nagligibla. Wyrtki (1981) haa axaninad tha avallabla avldanoa and 
concludad that upwalling valocitiaa at tha aquator ara of ordar 10*^ ai a~^ . 

Tha aquilibrivon aolution diacuaaad briafly abova with a • 3 , 6 1 

corraaponda roughly to thla cast and rasults in |6^ > d) <■ 3 K , axactly tha 
valua usad by Wyrtki in calculating tha nat horizontal haat tranaport away from 
tha equatorial "cool tongue". It ahould be pointed out, however, that horiaontal 
cold water advection is also inportant at the equator (as noted by Nyrtki) , 
contrary to the assumptions of the simple equilibrivnn model. 

The global heat balance for an area A of the cpol tongue is, retaining 
the entrainment term in Eq. (62) with definitions of and as before. 


— 2^ dxdy-(X hSu'nds- 
^w^w J 


■ - V - - V 


According to the estimates of Wyrtki (1981), the total volume tranaport 

transporting away much of the retained heat in the surface layer is of 
6 3-1 

order 50 x 10 m s over a 10,000 km long piece of the cool tongue in the 

equatorial Pacific. The corresponding kinematic heat transport, with 

(6, - 9. ) ■ 3 K is 150 X 10* K s ^ , or a large fraction of the world- 


wide poleward transport of heat by the ocean. 
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On* *sp*ct of *qu*torl*l h**t r*t*ntion has b**n Insufflciontly enphMlzod 

6 3-1 

in th* lltoratur*: th* 50 x 10 n a of h*at*d watar moving polaward 
constitutes a warm water mass formed in the equatorial cool tongue. The 
equivalent amount of water transport moving toward the equator to svqpply th* 
upwelling %iater has an average teng>eratur* scxn* 3 K lower, i.e. it is a 
different water mass, not a recirculating one as is sonetimes tacitly implied 
in discussions of the meridional circulation near the equator. At about 3*N 

6 

and S there is a total poleward moving water mass of order M ■ x 50 x 10 
and a similar equatorward moving mass, with a teiigmratur* some 3 K lower, all 
within 300 m or so of the surface. The poleward moving surface layer is 
strongly convergent beyond 3 *n or S because the Ekman transport drops sharply. 
How this warm water mass is incorporated into the equatorial near-surface 
circulation, and by what mechanism its replacement is supplied to the v^elling 
region, are two questions to which we have no answers at all at present. 
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Fig. 1 


Pig. 2 


Fig. 3 


Pig, 4 


Pigurg loggndi 

SchMiB of globftl oeoanic hMt transports, "cold” and "vam* 
watar mass formation in tha polar and aquatorial ocaans 

Intaracting atmospharic and ocaanic mixad layars and variablas 
of intarast in tha tharmodynamic modal 

Typical profilas of tamparatura and htanidity in tha atsiospharic 
mixad layar and abova, from Bunkar at al. (1949). Wind spaad 
was 7.5 m s saa minus air tamparatura diffaranca 0.5 K 

Annual cycla of surfaca haat fluxas and tamparatxiras : 

H nat radiant haat gain by tha ocaan, H sansibla, 

H, latant haat transfar to tha air, H • - H. 

L n r 8 L 

nat ocaanic haat ratantion, t^ saa surface, t^ air, dry 
bulb, t. wat bulb tamparatura. From Bunkar (1976), based 

d 

on more than 12,000 Observations in the trade wind region 
(23*N 52*W) 
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